Chemical equilibration of U(l) charged particles in a thermal environment 
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We discuss chemical equilibration of particles carrying non-vanishing quantum numbers related 
with U(l) internal symmetry. We construct the transport equation for the time evolution of particle 
multiplicities and their probability functions. The solution of these equations is obtained in different 
limiting cases. It is argued that a U(l) charged particles, dependent on thermal conditions inside a 
' ■ fireball, approaches different equilibrium limits. The differences between kinetics of abundantly and 

, rarely produced particles are explained. 
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. , I. INTRODUCTION 

> . 

' One of the crucial questions being addressed in the context of heavy ion collisions is the equilibration of the high 
density QCD medium created during these collisions jlj-Q. This question has been considered on different levels: by 
analysing physical conditions required for a perturbative QCD medium to reach equilibrium ||] or by studying the 
level of particle equilibration in the final state j^J^. From the theoretical point of view to discuss equilibration one 
fSJ ' needs to formulate the kinetic equation for particle production and evolution. In the partonic medium this requires, 
^ in general, the formulation of a transport equation involving colour degrees of freedom and a non-abelian structure of 
QCD dynamics In the hadronic medium, on the other hand, one needs to account for the charge conservations 

|, related with U(l) internal symmetry [pd|-p^. 
I ' In this article we will discuss the formulation of kinetic equations for U(l) charged particles produced in a thermal 
^-H [ environment. We will indicate the importance of the conservation laws in the time evolution and chemical equilibration 
T-H 1 of particle multiplicities and their probability distributions. In particular we will argue that the constraints imposed 
' by U(l) charge conservation are of crucial importance for rarely produced charged particle species. 

^ ■ 

^ . II. KINETIC MASTER EQUATION FOR PROBABILITIES WITH U(l) SYMMETRY CONSTRAINTS 

• To study equilibration in a hadronic medium one introduces a kinetic model that takes into account the production 
^ I and annihilation of particle-antiparticle pairs cc carrying U(l) quantum numbers. It is also assumed that particles c 
are produced according to a binary process ab cc and that all particle momentum distributions are thermal and 
described by Boltzmann statistics. The U(l) charge neutral particles a and b are constituents of a thermal fireball of 
' temperature T and volume V. We will consider the time evolution and equilibration of particles c inside a thermal 
fireball, taking into account the constraints imposed by U(l) symmetry. First we formulate a general master equation 
for the probability distribution of particle multiplicity in the medium with vanishing net charge and consider its 
properties and solutions. We will then discuss two limiting cases of abundant and rare particle production. Finally 
the rate equation will be extended to the situation where there is a net U(l) charge in the thermal fireball. 



- 1 — \ 

X 



A. Rate equation for particle multiplicity distribution 



Consider PN^{t) as the probability to find Nc particles c, where < iVc < oo. This probability will obviously change 
in time owing to production ab — *■ cc and absorption cc — *■ ab processes. The probability P/v^ tends to increase in 
time, following the transition from Nc — 1 and Nc + 1 states to the Nc state. It also tends to decrease since the state 
Nc makes transitions to Nc + 1 and Nc — 1 (see Fig. 1). The transition probability per unit time from Nc + I Nc is 
given by the product of the probability L/V that the single reaction cc ab takes place multiplied by the number of 
possible reactions which is, {Nc + '^){Nc + !)• In the case when the U(l) charge carried by particles c and c is exactly 
and locally conserved, that is if Nc + Nc = 0; then this number is just {Nc + 1)^. Similarly the transition probability 
from Nc ^ Nc + 1 is described by G{Na) (Ni,) /V, where one assumes that particles a and b are not correlated and their 
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multiplicity is governed by the thermal averages. One also assumes that the multiplicity of a and h is not affected by 
the ab — s- cc process. 
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FIG. 1. Schematic view of the master equation for the probability Pn{t) due to ab — » cc and the inverse process. 
The master equation for the time evolution of the probability P/v^ (''") can be written in the following form pd| : 



dT 



G 



L 



V 



V 



V 



V 



(1) 



The first two terms in Eq. (|l|) describe the increase of Pn^ (t) due to the transition from Nc — 1 and Nc + 1 to the Nc 
state. The last two terms, on the other hand, represent the decrease of the probability function due to the transition 
from Nc to A^c + 1 and Nc — 1 states respectively. 

Multiplying the above equation by Nc and summing over Nc, one obtains the general kinetic equation for the time 
evolution of the average number (Nc) = J2n =o ^cPN^iT) of particles c in a system. This equation reads: 



(2) 



The above equation cannot be obviously solved analytically as it connects particle multiplicity (Nc) with its second 
moments (Nc). However, the solution can be obtained in two limiting situations: i) for abundant production of c 
particles, that is when (Nc) 3> 1 or ii) in the opposite limit of rare particle production corresponding to (Nc) <C 1. 
Indeed, since 



(n!) = {Ncf + {sn!), 



(3) 



where (SN^) represents the fluctuations of the number of particles c, one can make the following approximations: 
i) for (Nc) » 1 one has, (Nc) ~ (Nc)^, and Eq. (2) obviously reduces to the well known form 



however, for rare production, particles c and c are strongly correlated; then, 
ii) (Nc) < 1 and (Nc) ~ (Nc); consequently Eq. (2) takes the form: 



(4) 



(5) 



where the absorption term depends only linearly, instead of quadratically, on the particle multiplicity. 

From the above discussion it is thus clear that depending on the thermal conditions of the system (that is its volume 
and the value of the temperature) we are getting different results for the equilibrium solution and the time evolution 
of the number of produced particles c. This is very transparent when solving the rate equations (4) and (5) . 

In the limit when (Nc) ^ 1, the standard Eq. ^ is valid and has the well known solution: 

(iV,)(r) = (7V,),qtanh(r/ro), (6) 
where the equilibrium value (-/Vc)cq of the number of particles c and the relaxation time constant Tq are given by: 
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(^c)cq = V^ , To = -^, (7) 

respectively, with e = G{Na){Na)/L. 

In the particular case when the particle momentum distribution is thermal, the gain (G) and loss (L) terms are 
just the thermal averages of the absorption and production cross sections with 

G ^ (jVc)cq(A^s)oq . . 

L (iVa)cq(iVfc)eq' 

where we have employed the detailed balance relation between production aab and absorbtion CTcc cross section for 
ab cc processes. 

Assuming a Boltzmann particle momentum distribution, the equilibrium average number of particles k is obtained 
after momentum integration as: 

(^fc)eq = ^VTmlK2{mk/T), (9) 

where d/c's denote spin-isospin degeneracy factors, rrifc the particle mass, and K2 is the modified Bessel function. 

This is a well known result for the average number of particles in the Grand Canonical (GC) ensemble with respect 
to U(l) charge conservation. The chemical potential, which is usually present in the GC ensemble, vanishes in our 
case, because of the requirement of U(l) charge neutrality of a system. Thus, the solution results in the expected 
value for the equilibrium limit in GC formalism where the charge is conserved on the average. 

In the opposite limit, where {N^) ^ 1, the time evolution of particle abundance is described by Eq. (5), which has 
the following solution: 

(iV,)C(r) = (A^,)C (l-e-/-o^), (10) 

with the equilibrium value and relaxation time given by 

{N,)^^ = e, = ^. (11) 

The above result is the asymptotic limit of the canonical C formulation of conservation laws ||l^. Here the charge 
related with U(l) symmetry is exactly and locally conserved, contrary to the GC formulation where this conservation 
is only valid on the average. 

Comparing Eq. (Q) with Eq. (|ll|), we first find that, for {Nc) ^ 1, the equilibrium value in the canonical formulation 
is far smaller than what is expected from the grand canonical result as 

Secondly, we can conclude that the relaxation time for a canonical system is far shorter than what is expected from 
the grand canonical value as 

= To(^c)eq « To, (13) 

since in the limit (ii) the equilibrium value {Nc)^^ <C 1. 

The above discussion shows the importance of the canonical description of quantum number conservation if the 
U(l) charged particle multiplicity is small. We also see that the volume dependence differs in the two cases. The 
particle density in the GC limit is ^-independent whereas, in the opposite case, the density scales linearly with V. 

We note that (i) and (ii) limits are essentially determined by the size of {6N^), the fluctuation of the number of 
particles c. The grand canonical results correspond to small fluctuations, i.e. (SN^) / (Nc)'^ <C 1, while the canonical 
description is relevant if the fluctuations are large, i.e. {6N^) / (Nc)'^ ^ 1. 



III. MASTER EQUATION - GENERAL SOLUTION 

In the previous section, we formulated a general master equation (1) for the probability to find the number of 
U(l) charged particles being produced by the binary ab — > cc process. The approximate solutions of Eq. (1) were 
considered in two limiting cases. For a large number of produced particles carrying U(l) charge, the time evolution 
and the equilibrium limit of Eq. (1) correspond to the GC result. We have also indicated the difference between 
GC and asymptotic C results on the level of rate equations for particle multiplicity. This difference is particularly 
transparent when comparing master equations for probabilities. 
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A. Master equation for abundant U(l) charged particle production 



For abundantly produced particles c and c through ab — > cc, process we do not need to worry about strong particle 
correlations due to U(l) charge conservation. This also means that instead of imposing U(l) charge neutrality 
conditions through Nc — Nc = 0, one assumes conservation on the average, that is (Nc) — (Nc) = 0. In this case the 
master equation (1) can be simplified. 

In the derivation of Eq. (1) the absorption terms proportional to L were obtained by constraining the U(l) charge 
conservation to be local and exact. For the conservation on the average, the transition probability from Nc to the 
(Nc — l) state is no longer proportional to {L/V)Nc but rather to {L/V)Nc{Nc), since the exact conservation condition 
Nc — Nc is no longer valid and the number of c particles can only be determined by its average value. In the GC 
limit, the master equation for the time evolution of the probability fW^(T) takes the following form: 

= ^{Na){Nb)PN^-l + ^ANc+ 1){Nc)Pn^ + i 
(It V V 

-^{Na){Nt,)PN^ -^Nc{Nc)Pn^. (14) 

Multiplying the above equation by A^c) summing over Nc and using the condition that (Nc) = {Nc), one recovers 
Eq. (4), the rate equation for {Nc) in the GC ensemble. The above equation is thus indeed the general master equation 
for the probability function in the GC limit. Comparing this equation with the more general Eq. (1), one can see that 
the main difference is contained in the absorption terms, which are linear instead of quadratic in particle number. 

Equation (14) can be solved exactly. Indeed, introducing the generating function g(x,T) for Pn^i 

oo 

g(:r,r)= ^ x^=P^,(t), (15) 

the iterative equation (14) for probability can be converted into a differential equation for the generating function: 

^^ = |Vi(l-a.)[ff'-Vi5], (16) 

with the general solution [^ : 

g{x, t) = go(l - xe-^)e^<'--)(^-'-^\ (17) 

where g' = dg/dx, f = {Ly/e/V)T and y/e = {Nc)eq given by Eq. (9). 

One can readily find the equilibrium solution to the above equation. Taking the limit r = oo in Eq. (17) leads to 

5cq(x)=e-^(i--), (^8^ 

with the corresponding equilibrium multiplicity distribution: 

This is a Poisson distribution with averaged multiplicity ^/e. 



B. Equilibrium solution - general case 

The master equation (14), describing the evolution of the probability function in the GC limit, could be solved 
exactly. The general equation (1), however, because the quadratic dependence of the absorption term, can only be 
solved numerically. Nevertheless the equilibrium result for particle multiplicity can be given. 

Converting Eq. (1) for Path's into a partial differential equation for the generating function 

oo 

g(x,r)= ^ :r^=P^^(r). (20) 
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one finds jTll 



— — ^y{^'x){x9 +9 -<^9)- (21) 

The equilibrium solution gcq{x) thus obeys the following equation: 

a^5eq + .9cq - e.9cq - 0. (22) 
By a substitution of variables {x — y^/4/e), this equation is reduced to the Bessel equation, with the following solution: 

1 



cq(-^) = , ^, /o(27ii), (23) 



where the normalization is fixed by g{l) = ^ Pn^ = 1- 

The equilibrium value for the probability function P^^ is now expressed from Eqs. (20-23) as: 

^^-'^'^^ 7o(VWF' ^^^^ 

We note that the equilibrium distribution of the particle multiplicity in not a Poissonian. This fact was first 
indicated in equilibrium studies in . In our case this is a direct consequence of the quadratic dependence on the 
multiplicity in the loss terms of the master equation (1). A Poisson distribution is obtained from Eq. (24) if ^/€ ^ 1, 
that is for large particle multiplicity. In this case the C ensemble coincides with the GC approximation. 

The result for the equilibrium average number of particles c can be obtained as: 

The equilibrium solution presented above coincides with the expected result for the particle multiplicity in the 
canonical ensemble with respect to U(l) charge conservation Thus, the rate equation formulated in Eq. (1) is 
valid for any arbitrary value of (Nc) and obviously reproduces the standard grand canonical result for large (Nc). 
We can study the chemical equilibration of U(l) charged particles following Eq. (1) independently from the thermal 
conditions in a system. 



C. Master equation for net U(l) charge in thermal environment 

So far, constructing the evolution equation for probabilities, we have assumed that there is no net U(l) charge 
in a system. In the application of the statistical approach to particle production in heavy ion collisions, the above 
assumption is only justified when the initial state is U(l) charge neutral and when considering particle yields in 
full phase space. However, because of experimental limitations, one often deals with data in restricted kinematical 
windows. Here the overall U(l) charge is no longer zero and a generalization of the above master equation is required. 

In the following we construct the evolution equation for P^ (t) in a thermal medium assuming that its net U(l) 
charge S is non-vanishing. The presence of the net charge requires a modification of absorption terms in Eq. (1). 
The transition probability per unit time from the Nc to the Nc — I state was proportional to {L/V)NcNc. Given an 
over all net charge S the exact U(l) charge conservation implies that Nc — Nc — S. The transition probability from 
Nc to Nc — 1 due to pair annihilation is thus {L/V)Nc{Nc — S). Following the same procedure as in Eq. (1) one can 
formulate the following master equation for the probability P^ (t) to find Nc particles c in a thermal medium with a 
net charge S: 

= -{Na){m)Pl_, + -{Nc + l)iNc + 1 - S)Pl^, 

- ^{Na){m)Pl - ^Nc{Nc - S)Pl, (26) 

which obviously reduces to Eq. (|^) for 5 = 0. 

To get the equilibrium solution for the probability and multiplicity, we again convert the above equation to the 
differential form for the generating function g^{x,T) = X)?^^=o ^^"-^Wc (''')■ 
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^^^^ = ^{l-x) {xg[^ +g'sil-S)~ egs) . (27) 
In equilibrium, dg^{x, T)dT — and the solution for g^^ can be found as follows: 

S /2 

9c<iix) = j^^^Is{2V^), (28) 

where the normahzation is fixed by g{l) — '^Pn — 1- 

The master equation for the probability to find Nc antiparticles c, its corresponding differential form and the 
equilibrium solution for the generating function can be obtained by replacing S* — > — S* in Eqs. (26-28) 

The result for the equilibrium average number of particles (Nc) cq and antiparticles {Nc)cq is obtained from the 
generating function using the relation: {Nc)cq = g'{^)- The final expressions read: 



The strangeness conservation is explicitly seen by taking the difference of these equations, which yields, the net value 
of the U(l) charge S. 

The results of Eq. ( p9| ) were previously derived from the equilibrium partition function by using the projection 
method The master equation derived here allows a study of the time evolution of particle production in a 

thermal medium and the approach towards chemical equilibrium. This equation explicitly accounts for U(l) charge 
conservation for both strongly correlated and uncorrelated particles. 

The results presented here can be extended to the more general case where there are different particle species 
carrying the U(l) charge being produced in a thermal environment p2|. 



IV. CONCLUDING REMARKS 



We have discussed the kinetics and chemical equilibration of particles carrying quantum numbers related with U(l) 
symmetry. Our analysis was restricted to the particular case of one kind of particle and antiparticle being produced 
by the binary process. However, this example is of physics interest as it can be applied to the study of such problems 
as the equilibration of heavy quarks in QCD plasma due to gluon-gluon fusion or strangeness and baryon production 
due to baryon-baryon and meson-meson interactions. The last process is of relevance in the application of statistical 
models to the description of particle production in low energy heavy ion collisions . We have discussed 

the procedure that allows the construction of the rate for time evolution of particle multiplicity and its probability 
distribution. Our result is quite general and applicable, as well as for abundantly and rarely produced U(l) charged 
particles. We have indicated the important differences between two limiting thermodynamical situations. In the 
first case the system approaches the equilibrium limit corresponding to the grand canonical approximation, in the 
second that corresponding to the canonical exact value. Of particular interest are the master equations for particle 
multiplicity distributions, which can be used to study not only the time evolution of the particle number but also 
all their higher moments This allows, for instance, a study of the time evolution of U(l) quantum number 

fiuctuations ||l^ , |l3|| , which has been proposed as a possible signal for QCD phase transition fl^ . 
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